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We present the results of a search for a very light CP-odd Higgs boson a? originating from 
top quark decays t — s> H ± b — > W ±{ 'alb, and subsequently decaying into t + t~. Using a data 
sample corresponding to an integrated luminosity of 2.7 fb" 1 collected by the CDF II detector in 
pp collisions at 1.96 TeV, we perform a search for events containing a lepton, three or more jets, 
and an additional isolated track with transverse momentum in the range 3 to 20 GeV/c. Observed 
events are consistent with background sources, and 95% C.L. limits are set on the branching ratio 
of t — > H ± b for various masses of H and a?. 

PACS numbers: 12.60.Fr, 12.60. Jv, 14.65Ha, 14.80.Da 
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The Higgs boson is the last unobserved particle of the 
standard model (SM) [T]. In the SM the Higgs boson 
mass is unstable to quantum corrections. This problem 
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is naturally solved in supersymmetric models [2 . In these 
theories the Higgs boson sector is more complicated. The 
minimal supersymmetric extension of the standard model 
(MSSM) contains five Higgs bosons: a light and a heavy 
CP-even Higgs (h and H), a CP-odd Higgs (A), and 
a pair of charged bosons (if )■ The next-to- minimal 
supersymmetric model (NMSSM) 3 further extends the 
MSSM to include an additional CP-even and CP-odd 
neutral Higgs bosons. In the NMSSM the lightest CP- 
odd Higgs can be below the bb threshold, so that the 
a® boson decays only into t + t~, cc, or gg. 

The existence of the very light boson has two im- 
portant implications. First, the decay mode of the SM- 
like Higgs h — > a®af{ becomes dominant and other SM 
decay rates are decreased, so that the SM-like Higgs bo- 
son avoids the LEP II direct limit g] . The light SM-like 
Higgs helps to solve the naturalness and fine-tuning prob- 
lems arising in the MSSM [3], In addition, the charged 
Higgs boson must not be much heavier than the W bo- 
son, which helps to reconcile apparent discrepancies in 
the LEP lepton universality measurements [S]. Such a 
charged Higgs boson could appear in top quark decays 
t — > H ± b, escaping current limits [3] due to a new open 
decay mode H ± — > W ± a®, which has not been inves- 
tigated before. This motivates a search for bosons in 
decays of top quarks. 

In the pp collisions at the Fermilab Tevatron the top 
quarks are produced mainly in pairs, and within the SM 
almost always decay into a W boson and a 6-quark. The 
NMSSM scenario considered above differs from the SM 
process by the presence of one or two a° bosons in the fi- 
nal state. As the a® boson decay products are expected to 
have low momenta, these could remain undetected with- 
out affecting the measurements of the ti cross section and 
properties of the top quark. 

In this Letter we report on the first search for a light 
CP-odd Higgs boson a? in decays of top quarks through 
the intermediate charged Higgs boson t —> H ± b —> 
W ±{ ' a\b assuming a\ — > t + t~. We analyze a data set 
corresponding to an integrated luminosity of 2.7 fb _1 of 
pp collisions at y/s = 1.96 TeV collected by the Collider 
Detector at Fermilab (CDF II) [7], searching in candidate 
ti events for the presence of \ow-pt tracks [5] that could 
be attributed to r-decay products. 

We select candidate it events using criteria developed 
for a it cross section measurement [9]. The data events 
used in the analysis are collected by triggers that iden- 
tify at least one high-p^ electron or muon candidate us- 
ing the online data acquisition system. Subsequently, 
each event is required to have a single isolated e or \i 
withp T > 20 GeV/c and |r?| < 2.0(1.0) for e(/x) USE]. 
We require missing transverse energy fir > 20 GeV [12], 
as evidence of a neutrino from the W-boson decay, and 
at least three jets with Er > 20 GeV and |7y| < 2.0, 
reconstructed using a fixed cone algorithm of radius 



R = y/{Arj) 2 + (A^) 2 = 0.4 [13 . Backgrounds to tt pro- 
duction are reduced by requiring at least one of the jets 
to be identified as a 6-quark candidate using the presence 
of a displaced secondary vertex [14] , and by requiring the 
scalar sum of the transverse energy of the lepton, fir, and 
jets (Hr) to be above 250 GeV. We observe 1052 events 
passing these selection criteria, which define a pre-signal 
sample. 

The main contribution to the selected sample of 
events comes from ti production, which we model us- 
ing the pythia 6.216 Monte Carlo (MC) generator [T5] 
for both SM and new physics top quark decays, assuming 
m t = 172.5 GeV/c 2 . We use an ALPGEN 2.13 [16] matrix- 
element generator interfaced to pythia 6.325 for model- 
ing W+jets and Z/ r y*+ jets production. Other sources 
of events in the pre-signal sample include diboson pro- 
duction (WW, WZ, ZZ) modeled with pythia 6.216 MC 
generator, and multi-jet QCD events modeled using a 
data-driven approach described in [T7] . The detector re- 
sponse in all MC samples is modeled by a GEANT3-based 
detector simulation [T5] , 

We search for r-leptons from a\ boson decays in ti 
candidate events by looking for at least one low-pr (3 
GeV/c < pr < 20 GeV/c) track in the central detector 
region \rj\ < 1.1. The track must be well-measured, i.e. 
it should have a sufficient number of hits in the track- 
ing chamber. To ensure that the track is consistent with 
being produced in a pp collision, the distance of closest 
approach of the track with respect to the beam axis is 
required to be small. The track must also originate from 
the same pp interaction as the isolated lepton by requiring 
that \ztrack — ziepton\ < 5 cm, where the z-coordinate cor- 
responds to the point of the closest approach to the nomi- 
nal beamline. To suppress backgrounds from jets the can- 
didate track is required to be isolated from other tracks 
in the event. We sum the pr of every well- measured track 
with pr > 0.5 GeV/c, including the candidate track, 
within a cone of AP < 0.4 around the candidate and 
with a z-position of origin within 5 cm of the candidate 
track z. We require that the ratio of the candidate track 
pr to the sum pr of tracks in the cone be at least 0.9. We 
also ensure that the track is not within AP < 0.4 of the 
lepton (e or fi) or a jet, used to define the ti candidate. 

The isolated tracks can arise from the hard parton- 
parton interaction producing the high-py lepton candi- 
date as well as from the "underlying event" (UE). In 
what follows, we include in our UE definition contri- 
butions from additional simultaneous proton-antiproton 
collisions. Non-UE isolated tracks come from physics 
processes where more than one lepton is produced but 
only one is identified, such as from Z/7*— > £ + £~ events 
where one lepton triggers the event, while the other one 
has a pr below 20 GeV/c, or is a r ± that leaves a low- 
Pr track. We use simulated events to model the track pr 
spectra corresponding to leptons from the vector boson 
decays. 
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We use data to model the characteristics of UE tracks. 
We analyze several different data samples to verify that 
the UE track pr spectrum is independent of the data 
source. We select Z boson candidates by requiring events 
to have two leptons ( "dilepton events" ) with an invariant 
mass consistent with a Z boson. We also study "lepton + 
jets" events by requiring only one lepton candidate, sig- 
nificant missing transverse energy, plus one or two jets. 
This data sample is dominated by events from W boson 
plus associated jets production. We also analyze several 
data samples of QCD multi-jet events collected by trig- 
gers that identify at least one jet. Each sample requires 
a different jet Et threshold. 

The fraction of events in which UE tracks satisfy our 
selection criteria is about 7.5%, and is consistent between 
samples within 15% relative uncertainty. The pt spectra 
of the isolated tracks for different data samples normal- 
ized to the same area are shown in Fig. [I] The track 
Pt spectrum for lepton+jets events is corrected by sub- 
tracting contributions from tracks corresponding to real 
leptons from Z/j*, diboson, and ti events. This is done 
by accounting for tracks originating from a or Z bo- 
son in our MC samples, where the reconstructed track 
is traced back to the charged particle in the decay chain 
of the vector boson. In Fig. [I] both corrected and uncor- 
rected track pt spectra are shown. After the correction 
the pt spectra agree with those from dilepton and QCD 
multi-jet events. 

Isolated Track p T - 1 Jet Events 



|5 Z_ Lepton + Jet 




| CDF Run II (2.7 fb ') 



1 U-l I I I I I I I I I I I I I I I I I I I LjJ 

5 10 15 20 25 

Track p T (GeV/c) 

FIG. 1: The isolated track pr spectra from lepton+jets, 
dilepton, and QCD multi-jet data samples for events with 
exactly one jet. 



We tested the data to determine whether there are any 
correlations between the pt spectra of isolated tracks and 
other parameters of the event. The only correlation we 
found was with the Ht of the event. We account for this 
correlation as described further in the text. A number of 
cross-checks that we performed include comparison of the 
UE track p T spectra for different QCD multi-jet samples, 
as well as a study of dependence on the jet multiplicity, 
the number of primary vertices, and presence of the b- 



tag in the events. We observed no statistically significant 
difference in the track pt spectra in these studies. 

We perform the search for t — s- H ± b — s- W ±( a\ b de- 
cays by fitting the observed isolated track p T distribution 
to the combination of the UE, non-UE SM, and the new 
physics signal track pt spectra. We use the UE track 
Pt distribution from QCD multi-jet data events to model 
the UE contribution, and allow the rate of UE tracks to 
float freely in the fit. For the MC- modeled background 
processes we consider isolated tracks only from the vec- 
tor boson decays. In case an event has more than one 
track satisfying the isolation and the track quality crite- 
ria, we select the track with the highest px- We use the 
UE track pt distribution measured in data to correct all 
MC track pt spectra to account for the probability of the 
highest-pT track to come from the underlying event. 

Prior to performing the fit in the signal region, we 
test our procedure in the control region defined by events 
with one lepton plus one or two jets. In this region the 
dominant non-UE contribution is from Z/7*— > e + e~ or 
Z/7*— > events, where the second lepton from 

Z/j* is not identified but passes our isolated track re- 
quirements, or Z/7*— > t + t~ events where one r decays 
leptonically and is identified as an electron or muon, and 
the other one is identified as an isolated track. The lep- 
ton track pt spectra from Z/j* decays are on average 
more energetic than the UE track pr, and assuming the 
UE-only hypothesis an excess of events is expected in the 
tails of the observed isolated track pt distribution. We 
test whether we are able to observe the excess of events 
attributed to Z/j*~^ £ + £~ events at the rate consistent 
with the expectation. The expected number of events 
from the Z/7*— !• £ + l~ process is obtained using the MC 
normalized to data under the Z mass peak. 

We perform a log-likelihood fit to the observed iso- 
lated track pt spectrum, with UE and Z/7* rates com- 
pletely unconstrained, and other MC-based contributions 
(top and dibosons) constrained to be within their theo- 
retical expectations. The fit is performed in the range 
3 < Pt < 20 GeV/c separately for events with one and 
two jets. The results of the fit are presented in Fig. [2] 
The extracted Z/j* contribution matches the expecta- 
tions within the statistical uncertainties. 

We then proceed to fit in the signal region, and employ 
the CL5 likelihood ratio test statistic [T§] to quantify the 
search results. The systematic uncertainties enter the 
CL5 fit as Gaussian-constrained nuisance parameters. 

The ti contribution is obtained from the data using 
the same technique as in the ti cross section measure- 
ment [9] . The uncertainty on the expected it event yield 
is due to the lepton identification and triggering (2%), 
6-tagging efficiency (5%), the jet energy scale (5%), the 
uncertainty in the estimate of backgrounds to ti (3%), 
and limited data statistics (6%) accounting for the total 
ti normalization uncertainty of 10%. 

The Z/^*+ heavy flavor contribution is normalized to 
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FIG. 2: The isolated track pr spectrum fitted for Z/7*+jets 
cross section in lepton + jets data events with one jet and 
two jets separately. In both cases the fit results are consistent 
with expected Z/'y* contribution. 



data under the Z mass peak, with the dominant uncer- 
tainty due to limited statistics of Z+ tagged jet events in 
data (8%). The uncertainty on the diboson {VV) back- 
ground is due to NLO calculations [20] and parton dis- 
tribution functions, taken conservatively to be 10%, lu- 
minosity (6%), and the jet energy scale (20%). 

Since we require the isolated track not to be within a 
reconstructed jet, the systematic uncertainty in the jet 
energy scale leads to events migrating to/from the sig- 
nal region, which results in an additional 3% uncertainty 
for all MC-based backgrounds. The uncertainty on the 
isolated track efficiency is 3%, and is determined using 
Z/j* events. 

The largest variations in the UE isolated track px spec- 
trum come from varying the Ht requirement for the can- 
didate sample. We use the shapes obtained from multi-jet 
data for very low and very high Ht, and interpolate these 
distributions to obtain an intermediate shape. The in- 
terpolation is parametrized with a Gaussian-constrained 
nuisance parameter and integrated into the fit. We allow 
the UE track pt distribution to change in the fit accord- 
ing to the value of this nuisance parameter |21) . 

The expected event yields in the signal region are pre- 
sented in Table[T] The first row in the table represents the 



805 215 11 19 - 1049 1052 
2.6±0.3 - O.liO.O 0.7±0.1 79±12 83±12 70 

TABLE I: Expected event yields in 2.7 fb -1 before the track 
requirement (first row), and with at least one isolated track 
(second row) with 3 GeV/c < pr < 20GeV/c. In the second 
row events are categorized based on the origin of the isolated 
track. The number of UE events is the expected number of 
events before the fit to the track pr spectrum. 



numbers of expected and observed events before the iso- 
lated track requirement. The second row shows the event 
yields after the isolated track requirement, where events 
are categorized based on the origin of the isolated track. 
The quoted event yield due to the UE corresponds to the 
expected rate, while the actual normalization is obtained 
from the fit to the isolated track px spectrum in data, as 
can be seen in Fig. [3] 

Figure [3] shows that the data are well described by SM 
background sources. We set 95% confidence level (C.L.) 
upper limits on the branching ratio of t — > H + b under the 
assumption that the branching ratios B{a^ — > t + t~) = 
100% and B(if ± ->■ W ± a\) = 100%. The expected and 
observed 95% C.L. limits as a function of to#± and m a o 
are shown in Fig. [3] For a given mass of CP-odd Higgs 
boson a? we exclude the branching ratios of B{t —> H + b) 
above the respective curve shown in the plot. For an 
a 5 boson with mass of 9 GeV/c 2 , we exclude a S(t -) 
H+b) > 0.20 at 95% C.L. for H+ masses between 90 and 
160 GeV/c 2 . These are the first limits on the branching 
ratio of t — > H + b in this decay mode. 

In conclusion, we have presented a search for non-SM 
top decays t H ± b W ±M a^b within the NMSSM 
scenario using a data sample corresponding to 2.7 fb -1 
of integrated luminosity in 1.96 TeV pp collisions. We 
see no evidence of r's from light Higgs a° decays, and 
set the world's first limits on the branching ratio of t —> 
H + b in this mode, assuming B(a\ — > t + t~) = 100% and 
B^H* -> W ± a° 1 ) = 100%. 
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